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Abstract 

We present a measurement of the cross section and the first measurement of the 
heavy flavor content of associated direct photon + muon events produced in hadronic 
colhsions. These measurements come from a sample of 1.8 TeV pp collisions recorded 
with the Collider Detector at Fermilab. Quantum chromodynamics (QCD) predicts 
that these events are primarily due to Compton scattering process eg — > C7, with the 
final state charm quark producing a muon. The cross section for events with a photon 
transverse momentum between 12 and 40 GeV/c is measured to be 46.8 ± 6.3 zt 7.5 
pb, which is two standard deviations below the most recent theoretical calculation. A 
significant fraction of the events in the sample contain a final-state bottom quark. The 
ratio of charm to bottom production is measured to be 2.4 it 1.2, in good agreement 
with QCD models. 
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Measurements of the inclusive spectrum of direct photons in hadron-hadron colhsions 
have provided important tests of quantum chromodynamics (QCD). Similar tests have been 
made with inclusive measurements of heavy flavor production (b and c quarks). The data 
and current perturbative QCD models do not agree well for both inclusive processes, giving 
insights into possible limitations of such models|I| [0]. Two previous measurements of the 
associated production of direct photons and charm quarks have provided checks of the charm 
quark content of the proton|§|[H through the Compton scattering process eg — >■ 07. We 
present here an analysis with an order-of-magnitude more events, collected by the Collider 
Detector at Fermilab (CDF), that provides a quantitative test of perturbative QCD. In 
addition, this new measurement is sensitive to the production of bottom quarks in association 
with the photon. 

The associated production of direct photons and heavy quarks in hadron collisions is 
expected to be a unique system for the study of the charm quark, with a 9:1 ratio of charm to 
bottom quarks in parton level QCD calculations]^] [Q. Typically in hadron collisions, heavy 
quarks are produced in the gluon-gluon initiated processes gg QQ and gg ^ gg ^ qQQ 
where one of the final state gluons splits into the heavy quark pair. In either case if the gluon 
energy is sufficiently larger than the bottom quark mass, the production of bottom pairs is 
approximately equal to that of charm pairs. In semi-leptonic decays of heavy quarks, the 
harder fragmentation function of the bottom quark leads to its dominance in these samples 
(for example a 1:4 ratio of charm to bottom in reference]^). The direct-photon Compton 
process, however, is proportional to the quark electric charge squared, which increases the 

1 



ratio of charm to bottom by a factor of 4. In addition, the intrinsic bottom quark content 
in the proton is 60% smaller than the charm quark content in our kinematic region. The 
combination of the quark charge coupling and the different proton content means the charm 
quark is expected to play the larger role in direct photon events. In this Letter we present the 
first measurement of the charm and bottom composition of direct photon events in hadronic 
interactions. 

The data for this analysis are from an integrated luminosity of 86 pb^^ of pp collisions 
collected with CDF in the 1994-95 Tevatron collider run (Run lb). The CDF detector and 
its coordinate system has been described in detail elsewhere0]P|. The events in the photon 
data sample discussed in this paper triggered the experiment by satisfying the requirement 
of a photon and a muon candidate at the hardware trigger level, whereas in the previous 
measurements only the photon candidate was required by the trigger. This allowed a lower 
transverse momentum Pt{= Psin{6)) threshold, in this case 10 GeV. A photon candidate is 
selected by requiring a cluster of energy in the central electromagnetic calorimeter \ri^\ < 0.9, 
with no charged tracks pointing to the cluster. The clusters are required to have a photon 
Pt between 12 and 40 GeV and to be isolated, with less than 1 GeV of additional transverse 



energy in a cone of AR = y^A(f^~+~Arf = 0.4 around the cluster. Additional photon 
cuts were used which were identical to those used in the Run la CDF inclusive photon 
analysis|TO|. Muon candidates were selected by requiring a match between a charged track 



with Pt > 4: GeV/c in the central tracking chamber and a track in the appropriate muon 
system. For \ri^\ < 0.6 muon candidates are required to be identified in both the CMU and 



in the Central Muon Upgrade system (CMP), which is behind an additional 1-m thickness 
of steel. For 0.6 < l?]^! < 1.0 the muon candidate track was required to be reconstructed in 
the central muon extension system (CMX). All three muon systems are discussed in detail 
in reference [ffl||TT|. After the track matching requirement there are 3850 events with a direct 



photon candidate and a muon candidate. 

Photon backgrounds from 7r° and rj meson decays remain in the sample. They are 
subtracted on a statistical basis using the photon background subtraction "profile" method 



described in reference ||10|| . This method uses the transverse energy profile of the electro- 
magnetic shower as a discriminant between single direct photons and multiple-photon meson 
decays. The requirement of photon Pt < 40 GeV/c described above is necessary in order 
to use this technique. After subtracting these backgrounds, 1707±83 direct photons with 
a muon candidate remain. For comparison, the previous publication Q was based on 140 
events. 

Muon backgrounds from charged pion and kaon decays remain in this sample, as well 
as a smaller fraction of charged hadrons that do not interact significantly in the material in 
front of the muon detectors. These are estimated with the same technique as in the previous 
analysis. Starting with the parent inclusive photon + jet data sample, the 4-vector of each 
charged particle with Pt > 0.4 GeV/c is measured. Each track is passed into a detector 
simulation as a charged pion or kaon, with a n/K ratio of 60%/20%[|12|. The results of 
the simulation are passed through the muon reconstruction; the tracks passing all cuts form 
the sample used for the background estimate. Backgrounds from protons that penetrate 



the calorimeter are negligible. After statistically subtracting the muon backgrounds, we 
expect 724±89 direct photon events with a muon that is not from charged vr^'' or decay. 
These events are assumed to come from associated direct photon + heavy quark production 
with the heavy quark decaying into a muon. Figure 1 shows the number of signal and 
background events in five bins of photon P^". Note that the purity of the sample improves 
dramatically as the photon Vt increases, which is due to the improved rejection of neutral 
meson backgrounds. The measured purity is consistent with that in inclusive direct photon 



measurements ||T0 1 



The photon-muon cross section da"'^^ / dP^ is derived for these five bins in photon by 
dividing by the luminosity, 86 the photon P^ bin size in GeV/c, and the efficiencies for 
detecting the photon within \ri^\ < 0.9 and the muon within |?7^| < 1.0. These efficiencies 
include the detector acceptance within the relevant pseudorapidity range, but are defined 
after the photon or muon Pt cut. The efficiencies are measured by a combination of studies 
using Monte Carlo simulation and data[Q|TI]. The photon efficiency varies from 34% to 



38% with a small photon dependence, while the muon efficiency varies from 49% to 53% 
and depends slightly on the specific muon subsystem. The resulting photon-muon cross 
section is shown in Table 1 along with the statistical uncertainties. 

There are four significant systematic uncertainties on the direct photon + muon cross 
section: 1) 12% from the muon background subtraction, which mostly comes from the uncer- 
tainties in the estimated pion and kaon particle fractions; 2) 7% from the photon background 
subtraction uncertainty, estimated in the inclusive photon measurement; 3) 7% from the 



Photon Ft Bin (GeV/c) 


daf+f/dP^ (pb/GeV/c) 


PYTHIA (pb/GeV/c) 


NLO QCD (pb/GeV/c) 


12-14 


7.5±1.9 


3.4 


10.3 


14-17 


4.4±1.0 


2.3 


6.4 


17-20 


1.7±0.7 


1.5 


3.8 


20-26 


1.5±0.4 


0.9 


1.9 


26-40 


0.3±0.2 


0.3 


0.5 



Table 1: The measured photon- muon cross section and the predictions from PYTHIA and 
NLO QCD are tabulated in five bins of photon P^. 



uncertainty in the photon and muon cut efficiencies; and 4) 4.3% from the uncertainty in 
the CDF luminosity measurement, which is predominantly due to the uncertainty in the 
total pp cross section. These added in quadrature give an uncertainty that ranges from 16% 
to 20% as the photon P^ increases from 12 to 40 GeV. 

The photon-muon cross section is compared to two different QCD calculations of photon- 
muon production. The first calculation is that in the PYTHIA Monte Carlo, which only 
has the leading-order (LO) contributions to the photon+heavy quark cross section, but has 
the full parton shower and fragmentation effects. The CLEO heavy quark decay tables 
are used[0. The second calculation is a next-to-leading order (NLO) QCD photon+heavy 
flavor calculation which has additional processes not present at leading-order, of which 
gg ^ cc ^ ccq is the largest contributor, but which operates only at the parton level and thus 
does not include heavy quark fragmentation. For this we use the Peterson fragmentation 
model in PYTHIA and the CLEO decay tables. The probability of observing a 4 GeV 
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muon from this model of the decay is shown in Fig. 2 as a function of heavy quark Pt 
for both direct bottom and charm decays and for sequential decays of the bottom quark. 
The NLO QCD calculation uses the massless quark approximation, which is adequate for 
photon+charm since the scale of the process is well above the charm mass, but cannot 
be applied to the photon+bottom process. The photon+bottom part of this calculation is 
derived from a leading-order photon+bottom calculation including mass effects, to which a 
multiplicative factor ('K factor', K) is then applied to account for higher-order effects. We 
use K= 1.8, which is the calculated ratio of NLO/LO in inclusive bottom production in a 



kinematic range that is close to that of the current measurement ||14|| . The uncertainties in 
the NLO calculation have not been thoroughly studied, but it is likely the K factor is the 
largest component. If one assigns a 50% error to K, then the predicted total photon+heavy 
flavor NLO cross section has a 12% uncertainty. Variations of renormalization scale, parton 
distributions, and fragmentation functions all give variations of ~ 5%. 

The NLO QCD cross section, as well as the PYTHIA predictions, are compared to the 
data in Table 1 and Figure 3. The NLO cross sections are much larger than those calculated 
with PYTHIA, due to the inclusion of additional processes mentioned above. The shape of 
the data shown in Figure 3 matches both calculations, while the normalization is a factor 1.9 
larger than the PYTHIA prediction and a factor 1.45 smaller than predicted by NLO QCD. 
With a data normalization uncertainty of 16% and a 12% uncertainty in the NLO QCD 
calculation, the data lie about two standard deviations below the NLO QCD prediction. 
The shape of the photon+muon cross section, however, is better described by theory than 



is the inclusive photon Pt spectrum. In addition, measured inclusive heavy flavor cross 
sections have typically been a factor of 2 larger than NLO QCD predict ions we do not 
observe this in photon+heavy flavor production. 

The current data sample is large enough to study the ratio of charm quark to bottom 
quark production in association with the photon. A "jet" of charged particles is measured 
by using the muon candidate as a seed, and then clustering charged tracks in a cone of 
radius 0.7 in rj — cj) space around the muon. We use the transverse momentum of the muon 
with respect to the jet the variable to separate the charm and bottom fractions 

of the sample. The larger bottom quark mass leads to an enhanced P^K A maximum 
likelihood fit is performed to the distribution of P^?' from the photon + muon candidate 
sample to four template spectra: 1) photon+charm, 2) photon+bottom, 3) photon + fake 
muon, and 4) 7r° + X where X is a muon candidate from any source (charm, bottom, or 
fake muon). The first two templates are generated using the PYTHIA Monte Carlo [0. 
The model is checked by returning to the parent photon+jet sample and comparing the P^'' 
distribution of the highest P^ track in the jet to the same from PYTHIA. The modeling is 
quite good using default PYTHIA parameters. The third template, photon + fake muon, 
comes from the same data+simulation combination used to estimate the fake muons for the 
cross section. Its normalization in the fit is constrained with a Gaussian weight using the 
systematic uncertainty of 12% coming from the cross section measurement. The fourth 
template, 7r° + X, comes directly from the data as it is the component subtracted during 
the direct photon background subtraction. The systematics on the maximum likelihood fit 
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using these four templates come from uncertainties in the shape of each template. The sys- 
tematic uncertainty on the photon+charm quark and photon+bottom quark templates are 
estimated by varying the input PYTHIA parameters. The allowed range of the parameters 
was determined by comparing the P^?^ distribution in the parent photon+jet data sample 
with PYTHIA. This leads to an uncertainty of 0.2 in the charm/bottom ratio. The system- 
atic uncertainty on the two data-driven templates, photon + fake muon and 7r° + X, are 
dominated by the photon subtraction method uncertainty and is 0.15 in the charm/bottom 
ratio. The statistical uncertainties from the maximum likelihood fit are much larger than 
the systematic uncertainties, being 1.1 in the charm/bottom ratio. The result of the fit to 
the charm/bottom ratio of the sample is 2.4±1.2 (stat. + sys.), to be compared with the 
predictions of 2.9 by PYTHIA and 3.2 by NLO QCD. There is good agreement between 
data and theory in this ratio. The unique nature of this charm-enriched sample is con- 
firmed, as the measured charm/bottom ratio of 2.4 is much larger than the value of 0.2 
in the inclusive heavy fiavor samples in Reference[0]. The photon+muon P^f^ distribution 
after all events with either a fake photon or muon are subtracted is shown in Figure 4a, 
compared to the templates for photon+charm and photon+bottom. The templates are 
normalized to the data. The photon+muon data are clearly a combination of charm and 
bottom, with a distribution more like the softer charm distribution. For comparison the 
same distribution using the 7r° + muon events from the same sample is shown in Figure 4b. 
Although QCD predictions for this process do not exist, qualitatively one would expect a 
larger bottom content since this process would not have the matrix-element enhancement 
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due to the square of the electric charge of the quark, as does the photon+muon 'Compton' 
process. The data confirm this hypothesis, albeit with limited statistics. 

In summary, we present a measurement of direct photon plus associated muon production 
in hadronic interactions with an order-of-magnitude more events than the previous measure- 
ment. The measurement is an interesting combination of direct photon and heavy flavor 
physics, each of which has had difficulties in comparisons with NLO QCD calculations. The 
data agree in shape with the theoretical predictions, but fall below the theory in normaliza- 
tion by 2 standard deviations. The ratio of charm/bottom in the sample has been measured 
for the first time, and confirms the QCD expectation that the sample is very enriched in 
charm quarks compared to inclusive lepton samples in hadron collisions. 
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Figure 1: The number of signal and background events 
in each bin of photon Pt are shown for the j + /J, sample. 
In addition to the signal of a real photon and a muon 
from a heavy-quark decay, the three components of the 
background are shown: 1) fake photons plus real or fake 
muons, 2) real photons plus a muon coming from the 
decay of a charged pion or kaon, and 3) real photons plus 
a fake muon coming from the punchthrough of a charged 
pion or kaon. 
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Figure 2: The probability for producing a muon with 

Pt >4 GeV in the semi-leptonic decay of a heavy quark 

is shown as a function of the quark Pt- These were 

generated using the PYTHIA Monte Carlo program to 

simulate the Compton 7 + process. The number of 

muons produced in association with the photon, N(7 + /x 

+ X), includes the relevant branching ratios (shown in 

parentheses). The quark and muon pseudorapidity are 

required to be within \r]\ < 1.0. The contribution from 

direct and sequential decays of the bottom quark are 

shown, as well as the direct decay of the charm quark. 
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Figure 3: The measured 7 + cross section is shown as a 
function of photon Pt. The photon pseudorapidity is re- 
quired to be within \rj\ < 0.9 and the muon pseudorapid- 
ity is required to be within \r]\ < 1.0. There is an overall 
16% normalization uncertainty (not shown) on the data. 
The measurement is compared to QCD predictions from 
the PYTHIA Monte Carlo as well as NLO QCD calcula- 
tions. Both the total contribution from 7 -|- (6 -|- c) and 
the individual 7 -|- c contribution are plotted. 
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Figure 4: The Pf^^ method is used to measure the sample composition of 7+heavy quark 
(left figure), and 7r°+heavy quark (right figure). In both figures, the points are the data 
after subtracting the contributions from fake muons. The solid curves are the best fit to the 
data, the dashed curves are the expected P^^^ distribution of a sample that is 100% charm, 
and the dotted curves are the expected P^p^ distribution of a sample that is 100% bottom. 
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